Coulomb exchange and pairing contributions 
in nuclear Hartree— Fock— Bogoliubov 
calculations with the Gogny force 



M. Anguiano, J.L. Egido and L.M. Robledo 

Departamento de Fisica Teorica, Universidad Autdnoma de Madrid, 
E-28049 Madrid, Spain 



Abstract 

We present exact Hartree-Fock-Bogoliubov calculations with the finite range den- 
sity dependent Gogny force using a triaxial basis. For the first time, all contributions 
to the Pairing and Fock Fields arising from the Gogny and Coulomb interactions as 
well as the two-body correction of the kinetic energy have been calculated in this 
basis. We analyze the relevance of these terms in different regions of the periodic 
table at zero and high angular momentum. The validity of commonly used approxi- 
mations that neglect different terms in the variational equations is also checked. We 
find a decrease of the proton pairing energies mainly due to a Coulomb antipairing 
effect. 
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1 Introduction 



The mean field approximation is the backbone of many-body calculations 
in Nuclear Physics, either as a zero order approach to the problem or as a 
basis for more elaborated theories like the Random-Phase approximation or 
the Generator Coordinate Method. The modern and fast computer facilities 
have made possible the use of effective density dependent interactions, as the 
Skyrme [1] or the Gogny force [2] , in standard mean field calculations like the 
Hartree-Fock (HF) or the Hartree-Fock-Bogoliubov (HFB) approaches. 

In spite of the relative simplicity of the mean field approach, some additional 
approximations have been commonly used when effective forces were used 
in the mean field calculations in order to make the problem more tractable. 
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The complications in the calculations usually arise from the exchange terms, 
either from the Coulomb force or from some components of the nuclear force 
itself. Depending on the mean field approach (HF or HFB) and on the force 
(zero or finite range) different approximations have been used in the past. 
In the HF theory one has only one exchange term (the Fock term) which 
contributes to the HF field. Commonly, in this type of calculations, in order 
to retain the required simplicity of the model, the Fock term of the Coulomb 
force has been neglected, or treated in the Slater approximation. In ref. [3] the 
validity of different approaches for the Coulomb Fock contribution was tested 
in HF calculations finding the Slater approximation a good one. In the HFB 
approach, and for any two-body interaction, one has two exchange terms, the 
Fock term which contributes to the HF field and the particle-particle (pp) 
term which contributes to the pairing field. In the HFB approach one must 
distinguish between the Skyrme and the Gogny force. Since the Skyrme force 
is assumed to provide only the particle-hole (ph) part of the force, all terms 
of the pp type are neglected even those arising from the Coulomb force. With 
the Gogny force the situation is different because with its finite range it has 
been designed to provide both, the ph and the pp part of the nuclear force. 

A test similar to the one performed for HF in ref. [3] has never been done for 
HFB calculations. Exact HFB calculations with the Gogny force have been 
performed for spherical nuclei [2], but the effects of the approximation of 
neglecting the exchange terms never have been analyzed. For the more complex 
triaxial HFB calculations it has been assumed that the findings of ref. [3] for 
the HF approach, concerning the exchange terms, will also apply to this case. 
This point, however, never has been checked. 

HFB calculations with triaxial symmetry using density-dependent interactions 
like the Gogny force, used by different groups [4-6], are usually performed 
using the following approximations: 

• First, the contributions to the pairing field stemming from the center of 
mass correction. Coulomb and spin-orbit terms are neglected. This is moti- 
vated by the fact that their contribution, in the limit of axial deformations, 
is assumed (but not proved) to be negligible small [4]. Thus, under these 
conditions the pairing field arises only from the Brink-Boeker term. 

• Second, since the exchange (Fock) Coulomb term contribution to the HF 
field requires a large CPU time, its contribution to the energy is not taken 
into account or it is included in the Slater approximation. 

On the other hand in the HFB theory, pairing is an important degree of free- 
dom and it is important to analyze how the properties related to this part of 
the field change when some approximations are assumed. In this paper, for 
the first time, we calculate in the triaxial basis the contributions mentioned 
above . We shall investigate the effect of including all the terms in the fields 
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as compared with the calculations which are normally done for deformed nu- 
clei. Mainly the DIS parameterization of the Gogny force [8] is used in the 
calculations, though some calculations with the Dl [2,9] parameterization will 
be done. We shall study nuclei which have been analyzed by the approximate 
triaxial HFB, hke some nuclei in the rare earth region [5,10], in the superde- 
formed Aps 190 region [6] and in the actinide region [4]. In the same way, 
we shall also analyze these contributions in spherical nuclei in order to un- 
derstand how the results obtained depend on the basis used. A brief review 
of the method is given in Sect. 2. The effect of the different contributions in 
spherical nuclei is studied in Sect. 3.1 for the N= 50 region and in Sect. 3.2 
for the N= 82 y N= 126 regions. The rare earth region will be analyzed in the 
Sect. 4.1. The effect of the new terms in the nucleus ^^°Pu (actinide region) 
will be treated in Sect. 4.2. The behavior of a nucleus in the super deformed 
region A^i 190 will be studied in Sect. 5. Finally an study of energy surfaces 
with considerations of the exchange terms will be done in Sec. 6. A discussion 
and the conclusions will be presented at the end of the paper. 



2 Theory 



The HFB theory [11] unifies the self-consistent description of nuclear or- 
bitals, as given by the Hartree-Fock (HF) approach and the Bardeen-Cooper- 
Schrieffer (BCS) pairing theory into a single variational theory. As trial wave 
function an independent quasiparticle state 1$) is considered. This wave func- 
tion is a linear combination of independent multi-particle states representing 
various possibilities of occupying single-particle states. The quasiparticle op- 
erators are defined by 



with c[, Cfc the particle creation and annihilation operators in the harmonic os- 
cillator basis and U and V the Bogohubov wave functions to be determined by 
the Ritz variational principle. The corresponding equations have been solved 
by the Conjugate Gradient Method [12]. The Bogoliubov transformation (1) 
leads to wave functions 1$) that are not eigenstates of the particle number 
operator. Furthermore, to generate a wave function 1$/) corresponding to a 
rotational state with angular momentum I we shall use the cranking prescrip- 
tion. Therefore, in the minimization process one has to add Lagrange multi- 
pliers to satisfy the pertinent constraints. That means, the wave function 
is determined by the condition 

6{(!>i\H -cuX- XnN - XzZ\^i) ^0 (2) 
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with the Lagrange multipliers uj.Xn and Xz determined by the usual con- 
straints (4) = ^/(J+ 1), (iV) = N and (Z) = Z, in an obvious notation. 

In the second quantization formalism, the nuclear Hamiltonian is given by 

kik2 kik2k3k4 

where t is the kinetic energy and v the two-body interaction, v if antisym- 
metrized. The Wick theorem allows the evaluation of the expectation value of 
a two-body operator in a simple way: 

(4l42Cfc4Cfc3) = (4lCfc3)(42Cfc4) - (42Cfe3)(4lCfc4) + (4l 42 ) (Cfe4Cfc3 ) , (4) 

where we have introduced the shorthand notation = (6) for the 

expectation value of any operator O. The first term on the r.h.s. of eq. (4) is 
called the Hartree term, the second one the Fock term and the third one the 
pairing term. Using the expressions 

PklM = {Ck2'^ki_), /«fci,fc2 = (Cfc2Cfei), (5) 

for the density matrix p and the pairing tensor k, and 

Tk^ks = E Vk,k2k3k4Pk4k2, ^k,k2 = ^ E Vk,k2k3k4K,,,,, (6) 
A;2fe4 fe3fc4 

for the HF field F and the pairing field A, we can express the expectation 
value of the nuclear Hamiltonian by 

E={H) = Tr{tp) + hpr (Fp) - (A«;*) . (7) 

The second addend in this expression arises from the contributions of the 
Hartree and the Fock terms of eq. (4) while the third one stems from the 
pairing term of eq. (4). The two-body effective nuclear interaction v used 
in our calculations, see Appendix A, is made up from the Gogny force, the 
Coulomb (C) interaction and the two-body correction of the kinetic energy 
(TK). The Gogny force itself has the following terms: Brink-Boeker (BB), 
spin-orbit (SO) and density-dependent (DD) contributions. Then, the total 
HFB energy can be written as 

E = Tr{tp) + Vbb + Vso + Vdd + Vtk + Vc (8) 
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where we have written separately the different energy contributions. Each 
contribution, Vl, is calculated from the corresponding Hartree-Fock, F^, and 
pairing, A^, fields (L = BB, SO, DD, TK or C) and is given by 

Vl = \tt {VlP - Az,^*) = + V[ + V[ (9) 

where F^ and A^ are given by eq. (6) but considering, instead of the full 
interaction f , only the part L of the force. On the r.h.s. of this expression we 
have further separated Vl into the three contributions obtained from each of 
the terms of eq. (4), namely, the Hartree-, the Fock- and the pairing term. 

As mentioned above, in HFB calculations with the Gogny force and a triaxial 
basis several approximations have been done. In particular the contributions 
to the pairing field from the spin-orbit, two-body correction of the kinetic 
energy and Coulomb terms have been usually neglected in the self-consistent 
procedure. Furthermore, the exchange contribution of the Coulomb potential 
(Fock term) to the Hartree-Fock field has not been taken into account in most 
of the cases or treated in the Slater approximation. Thus, in the simplest HFB 
approximation, which we shall denote HFB^ (s for standard), the HFB energy 
used in the variational process is given by 

Es = Tr{tp) + Vbb + Vdd + ^Tr ( (F50 + Ttk + T^)p) (10) 

A first order perturbation theory correction to this approach, which we shall 
denote HFBs+, is to calculate the neglected terms just at the end of the itera- 
tive procedure and to add them to the total energy. We shall test, in different 
mass regions, how good this approach is with respect to the exact calculation 
(7), denoted HFBg (e for exact), and with respect to other approaches that we 
shall introduce later on. Details on the calculation of the terms neglected in the 
HFBs approximation but included in the exact one are given in Appendix B. 

The triaxial basis used in the calculations is spanned by harmonic oscillator 
states with quantum numbers {nx,ny,nz} which fulfill the condition, GxUx + 
ttyUy + ttzTiz < Nq. The coefficients a, are related to the axis of the nuclear 
matter distribution by = (qp)^^^, cty — q^/^p~'^/^ and — p^/^q~'^/^ [4]. 



3 Spherical nuclei 

In this section we concentrate mainly on the effect of the Coulomb contribu- 
tion in HFB calculations for spherical nuclei. We are interested to know the 
behavior of the pairing and total energies when different approximations for 
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the Coulomb terms are Tised. At the end of sect. 3.1, a short discussion is 
devoted to the contributions to the total energy of the pairing field of both, 
the spin-orbit term and the two-body correction of the kinetic energy term. 



3.1 N— 50 region 

In this subsection we shall study the medium heavy nuclei from the region 
= 50. We shall use the basis determined hj q = 1.0, p = 1.0 and Nq = 11.1. 
In the calculations the DIS parameterization of the Gogny force has been used. 
First we analyze in detail the Coulomb terms and ignore for the moment the 
SO and TK pairing exchange terms. 

The Coulomb interaction contributes to the energy with three terms : The 
Hartree (H / h) or direct term, the Fock (F / f) or exchange term and the pair- 
ing (P / p) or Bogohubov term. Each term can be considered either selfconsis- 
tently, i.e., taken into account in the variational process, or non-selfconsistently 
(perturbationally), i.e., ignored during the variation but its value being added 
to the total energy at the end of the variational process. In the first case we 
shall use capital letters (H, F or P) and in the second one small letters (h, f or 
p) to label the approximation used for the respective terms. For example the 
approximation Hfp means that the Hartree Coulomb term has been taken 
into account selfconsistently and the Fock and Pairing Coulomb terms only 
perturbationally. To understand which contributions are more important and 
the convenience of considering in the variational process the pertinent terms 
we shall investigate the following approximations: 

In the simplest calculation the three Coulomb terms mentioned above, i.e., 
even the Coulomb Hartree term, are neglected. In this approximation the 
energy, which we shall call reference energy, Eref, given by 



is considered in the minimization process. The next approximations will be 



We use squared brackets to indicate the terms considered perturbationally, the 
others, are used in the variational equations. Notice, therefore, that only the 




(11) 



Ehfp = Eref + [Vc] 
Eufp = Eref + + [Vc^ + V^] 

EhFp = Eref + V^ + VS + [VS] 
EhFP = Eref + Vc 



(12) 
(13) 
(14) 
(15) 
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latter ones determine the wave function. To remember the meaning of each 
approximation one just needs to look at the subindex of the energy and to 
keep in mind the convention introduced above about capital and small letters. 
We shall also investigate the Slater approximation used by several groups. In 
this case the energy looks hke 

EhFp = Eref + + VE + (16) 

In i^^f^p the tilde on F means that the Fock term has been calculated in the 
Slater approximation but it has been considered in the variational process. In 
these calculations the best one is, obviously, E.FP where all Coulomb terms 
are taken into account selfconsistently. 

In the upper part of Table 1 we present the binding energies of five spherical 
nuclei with N — calculated in the different approximations. The exper- 
imental binding energies are from ref. [13]. As expected from a variational 
method, the more terms are treated selfconsistently the deeper gets the en- 
ergy minimum. The largest energy decrease is provided by the inclusion of 
the Hartree term in the variational equations. The largest deviation of the 
E.FP energies from the experimental values is about 2 MeV, but already the 
selfconsistent treatment of the Hartree term provides a good approximation 
to it. We also find that the Slater approximation provides energies close to the 
Hjj) approximation. These effects can be more clearly seen in Fig. la where 
we plot the quantity SEapp, which represents the percentage of error of a given 
approximation with respect to the calculation HFP and is defined by 

5E^pp, = 100 X ^" V (17) 

rjHFP 

where Eapp. represents any of the approximations that we have considered. 
Here, we see that the hfb approximation, where none of the Coulomb terms 
is treated selfconsistently, deviates the most from the exact HFP results and 
that the deviation increases with the proton number. The other approaches 
provide a good approximation to the complete HFP. 

In the second section of Table 1 we show the proton pairing energies E^^p for 
the ground state of the same nuclei. The neutron pairing is zero since we have 
a major shell closure for these nuclei. These energies are defined by the expres- 
sion EPPp^ — — ^Tr (Ak*). One has to notice that from the three Coulomb 
terms only the pairing term gives a contribution to Epp^ , the Hartree and 
Fock terms only give an indirect contribution through the changes introduced 
in the wave function by their consideration in the variational equations. We 
have also included in the table E^^j-, which gives us the pairing energy with- 
out any Coulomb contribution. The difference between E^^j: and -E/^/p gives 
us an indication of the size of the pairing Coulomb energy, which is about 1 
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Table 1 

Binding energies, proton pairing energies, lowest two-quasiproton excitation ener- 
gies and two-protons separation energies (in [MeV]) for spherical N = 50 nuclei. 

MeV for all nuclei. When the direct Coulomb term is included in the HFB 
calculation, Hfp, we observe a slight decrease ( ranging from —0.150 MeV to 
—0.670 MeV) in the pairing energies with respect to the calculation hfp. The 
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inclusion of the Fock Coulomb term, HFp^ produces the opposite effect, we 
get a slight increase (up to 1.0 MeV) of the pairing energy, again with respect 
to hfp. Finally, treating all the contributions of the Coulomb force in a self- 
consistent way, HFP, we observe an important increase (in absolute value, a 
decrease) in the values of the proton pairing energy of up to 4 MeV for the nu- 
cleus ^°Zr. It is interesting to notice that though, for this nucleus, the pairing 
Coulomb term itself is about 650 keV large, when considered selfconsistently 
it grows up to 4 MeV. Concerning the Slater approximation, HFp, we find 
again that it provides results very close to the Hfp ones. We observe that the 
maximum effect of the three contributions occurs for Z — 40 { semi-magic nu- 
cleus). This can be easily seen in Fig. lb where a definition similar to eq. (17) 
has been used for the pairing energy. Comparing a) and b), we find a clear 
difference: the percentage of error in the pairing energy is much bigger than 
in the binding energy. We obtain that even the Hfp approximation gives us a 
good result for the total energy, while this approximation is not good for the 
pairing energy. As a matter of fact, an important increase ( > 1.5 MeV for 
all nuclei) is obtained for the pairing energy when we compare HFP with all 
other approaches considered. 

In the third section of Table 1 we present the lowest two-quasiproton energies 
calculated in different approximations. In this case we have also included the 
exact solution, HFBg. In line with the decrease of the pairing energies just 
discussed, we observe a decrease of the excitation energies as more terms are 
included in the variational procedure. This decrease is particularly strong for 
the pairing exchange Coulomb term. We may conclude from this discussion 
that neglecting the Fock and pairing Coulomb terms in the variational equa- 
tions does not affect significantly the total energy but it does infiuence the 
pairing properties of the wave function. Then, if we arc interested in proper- 
ties of spherical nuclei for which pairing correlations play an important role, 
like excited states or transition energies, one must be careful with the Coulomb 
exchange terms. Notice that not only the two-quasiparticle excited states will 
be affected but also any calculation based on a HFB solution, like the Ran- 
dom Phase Approximation, will also experience a lowering in the excitation 
energies of the low-lying collective states. We can see in Table 1 that the most 
important reduction in the pairing energy is due to the introduction of the 
Coulomb pairing term in a self-consistent way, though the contribution of this 
term itself is not very important, about 1 MeV for the considered nuclei. 

From the results just discussed one could conclude that the Coulomb inter- 
action produces an antipairing effect in nuclei. In order to understand this 
effect we shall now analyze in a closer detail the behavior of the most relevant 
Coulomb term concerning pairing correlations, namely the pairing Coulomb 
term. We have done an additional test for the nucleus ^°Zr. First of all, an 
exact HFB calculation (including the three Coulomb terms and the spin-orbit 
and two-body kinetic mass correction contributions to the pairing field, i.e.. 
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Fig. 1. Percentage of error in a) the binding energy and b) the pairing energy of 
various approximations for the Coulomb part related to exact calculations as a 
function of the proton number Z. 

HFBg) has been done for this nucleus. Then, we do a progressive switching off 
of the Coulomb pairing contribution by a parameter r]. For example, i] = 1.0 
means that we take into account the full Coulomb pairing selfconsistently, 
rj — 0.6 means that 60% of it is taken selfconsistently and 40% perturba- 
tionally and so on. In fig. 2a we display the proton-proton correlation energy 
versus the r] parameter for the nucleus ^^Zr. We find a linear behavior, which 
means that the change in the pairing energy due to the switching off' of the 
Coulomb pairing contribution is produced in a continuous way. The solution 
that we obtain for 77 = 1.0 is very similar to the solution for 77 = except for 
a decrease in the proton-proton pairing energy. 
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Fig. 2. a) Epp versus the r] parameter for the nucleus ' ^ Zr. b) Particle energies e in 
the canonical basis and c), d) occupation numbers i^{nlj) versus the rj parameter 
for the nucleus ^^Zr. 

In fig. 2b) we present the proton single-particle energies, s{nlj), outside the 
Z— 28 shell closure in the canonical basis, and in fig. 2c), d) the occupation 
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lylnlj) for each level (nlj), which is given by 

^("^^■) = ^TTT E(^l4mCn/,m|$) (18) 

In this figure we observe that the single-particle energies e barely change with 
rj. With respect to the occupation numbers, we observe that the Coulomb 
pairing contribution goes in the direction to fill up the levels O/5/2, IP3/2 and 
lpi/2 and to empty the level Og9/2- That means, the Coulomb pairing term 
tries to increase the occupancy of the levels below the Fermi level at the cost 
of decreasing the occupancy of the levels above. 

In the forth section of Table 1 we show the two-proton separation energies 
S^P{N, Z) = E{N, Z -2)- E{N, Z). As we can see in the table we obtain 
similar results in all approximations and the agreement with the experiment 
is good in all cases. We can understand this result if we consider that the 
binding energies are well described in all approximations and that that S"^^ is 
not an observable related with the pairing properties of the nuclei, which are 
really sensitive to the different approximations. 

Finally, to complete the study of the neglected terms we shall now investigate 
the contributions of the spin-orbit term and of the two-body correction of the 
kinetic energy to the pairing field. In the fifth section of Table 1 we show 
the binding energies Eso,tk obtained by adding these terms at the end of the 
minimization of eq. (15), i.e., 

Eso,tk = Ehfp + [Vio + V^^] (19) 



in a similar notation as introduced above. These terms are repulsive and a 
few hundred keV large. If these two terms are now taken into account selfcon- 
sistently we obtain the exact solution, i.e., HFBe, of the full HFB equations, 
eq. (8). The resulting energies are on the average 100 kcV smaller than if the 
terms are considered pcrturbationally. The last two entries of the table are the 
pairing energies -E'^o j^, defined as before. We observe that these terms provide 
an additional quenching of the pairing correlations. The behavior is similar to 
the one observed with the other terms, i.e., the largest quenching takes place 
for the semi-magical nucleus ^°Zr. For this particular nucleus we have calcu- 
lated the effect of taking into account selfconsistently only one term and the 
other one pcrturbationally, i.e., we have performed the following calculations 

Eso,tk — Ehfp + V^q + [V^fj^] (20) 

Eso,TK — Ehfp + Vtr + Wso]- (21) 

We obtain Efo,tk = -1-671 MeV, Eso,tk = -784.360 MeV and E^j^j^ = 
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— 1.355 McV, Eso,TK = —784.354 McV, here again capital letters mean self- 
consistent and small letters pcrturbational approaches. As we can see each 
term provides about half of the total effect. 



3.2 The 82 and 126 regions 



In this section we focus on the chains of the N= 82 and N= 126 isotones, again 
with neutron shell closure. Then, only proton pairing correlations are going 
to be important in this region. Two approximations are used in this case: the 
approximation HFBs, see eq. (10), which neglects all the terms mentioned in 
the introduction, and the exact calculation HFBg, see eq. (8), which takes 
into account all the contributions from the different terms of the force to 
the fields. In the HFB^ approximation, the missing terms are not added at 
the end of the variational procedure. In order to be able to compare the total 
binding energies in the HFEg approach with the exact ones one must add these 
terms perturbativaly, i.e, the HFBs_|_ approximation. In the N= 82 region, the 
basis is determined hj q = 1.0, p = 1.0 and A^o = H-l and in this case the 
parameterization Dl is used. In the N= 126 region, the basis q = 1.0, p = 1.0 
and A^o = 12.1 has been used and for the force parameterization we have 
chosen DIS. 

In the first two columns of Table 2 we display the total binding energies. 
We find that the approach HFB^.!. provides a very good approximation to 
the exact results. In the last three columns of the same table we show the 
proton-proton pairing energies obtained in the approximations HFE^, HFBs+ 
and the exact ones HFBg. Both approximations provide more proton pairing 
correlations that the exact calculation HFBg. We can see that, contrary to the 
total binding energy case, the perturbational calculation HFB^.!. is not a good 
approximation for the pairing energy. We do not find any noticeable difference 
between both mass regions or the Gogny parameterizations. 

From these results we may conclude that, for spherical nuclei, the approxi- 
mation of ignoring some terms of the hamiltonian, HFBs, presents some dif- 
ferences with the exact calculation. The approximation of considering these 
terms at least in first order perturbation theory, HFBs_|_, works well for observ- 
ables related to, or that can be obtained from, total binding energies. On the 
other hand, the wave function content may be different from the one obtained 
in the exact calculation. 
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Nucleus E,+ Ee E^i W 



s s+ e 



132 

50 Sn 


1 Ann OTK 


-1090.891 


-0.0000 


-0.0000 


-0.0000 




-1111.079 


-1111.236 


-7.4700 


-6.5881 


-4.7615 


54 ^6 


-1129.549 


-1129.784 


-11.7377 


-10.3805 


-7.4571 


138t3„ 
56 


-1146.233 


-1146.557 


-14.3629 


-12.7469 


-8.6469 


140 

58 


-1161.102 


-1161.511 


-15.7718 


-14.0409 


-8.8442 



206 XT 

80 


-1615.860 


-1616.038 


-6.6570 


-5.8068 


-4.2799 


208 pu 

82 


-1634.620 


-1634.639 


-0.0000 


-0.0000 


-0.0000 


210p„ 

84 


-1642.746 


-1642.885 


-6.1489 


-5.3940 


-3.9841 


86 


-1649.656 


-1649.878 


-10.5256 


-9.2758 


-6.6871 


214pa 
88 


-1655.273 


-1655.560 


-13.7754 


-12.1942 


-8.2135 



Table 2 

Binding energies and proton pairing energies (in [MeV]) for spherical N = 82 and 
N = 126 nuclei with the methods HFB, and HFBp for I = Oh. 



4 Normal deformed nuclei 



In this section we will study the effects neglecting some terms in the hamilto- 
nian for deformed nuclei. In particular we shall investigate their behavior at 
high angular momentum as well as the impact of the exchange terms on the 
neutron system. 



4-1 Rare earth nuclei 



In this region, we study six nuclei: the three Erbium isotopes I64,i66,i68gj.^ 
the three Ytterbium isotopes, i62,i64,i66y]3 As in the preceding subsection we 
shall investigate only the approximation HFBs (or HFB5+ when the last one 
does not make sense) to compare with the exact calculation HFBe and the 
experimental results when available. For these calculations we use the basis 
with q= 1.3, p= 1.0 and No = 11.1 and the DIS parameterization of the force. 

We shall first investigate the ground state properties of these nuclei. In the 
upper part of Table 3 we show the binding energies. We find that the ap- 
proximate calculation, HFB^.!., provides results very close to the exact ones. 
In the middle of the table we display the pairing energies calculated in both 
approaches. We find that the exact calculation provides proton pairing ener- 
gies much smaller, in absolute value, than the approximate calculation, HFBs. 
The global efi^ect that we observe in the HFBe calculations is a decrease of the 
absolute value of the proton pairing correlations compared to the approximate 
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1 fiA -M~, 


1 fifi.-M~, 


168 Er 


1 AO-K rt 

162Yb 


1 fiA T-t 

164 Yb 


166 Y|_, 




-1328.665 


-1343.424 


-1357.391 


-1305.810 


-1322.764 


-1339.332 


Ee 


-1329.176 


-1343.868 


-1357.777 


-1306.361 


-1323.425 


-1340.123 




-7.068 


-6.297 


-5.642 


-10.081 


-9.354 


-8.309 




-1.723 


-1.727 


-1.730 


-4.167 


-0.001 


-0.000 




-8.757 


-7.475 


-5.602 


-9.444 


-9.449 


-9.049 




-7.756 


-6.310 


-4.310 


-8.223 


-11.005 


-9.614 


Q. 


7.590 


7.781 


7.893 


6.073 


6.871 


7.563 


Qe 


7.552 


7.712 


7.836 


6.190 


7.602 


7.787 


Table 3 



Ground state binding energies, pairing energies (in [MeV]), and quadrupole mo- 
menta (in barns), and in the approximation (HFBs) and in the exact one (HFBe). 



calculation. For -"^^^Yb and -"^^^Yb these correlations even vanish. The behavior 
of the neutron pairing energy is different depending on the nucleus: for the 
last two nuclei {E^^'^l is increased by about 1 MeV and for the other nuclei 
it is decreased by approximately the same amount as compared to \EpP''^\. In 
general, we can say that the global effect of including the Coulomb exchange 
and pairing contributions in the self-consistent procedure is a strong reduc- 
tion of the absolute value of the proton pairing correlations in all the cases. 
The circumstance that the inclusion of the neglected terms influences more 
the proton pairing than the neutron one is obviously due to the Coulomb 
terms that do not affect the neutron system. With respect to the quadrupole 
moments we observe that the relative variation between both calculations is 
relative small, approximately 0.5% for all the nuclei except for ^^^Yb where Q 
increases in the exact calculation by about 10%. 

The fact that the total binding energies are very close in both approaches is 
the reason why the non self-consistent and perturbational approximation has 
been considered in the past as a good approximation in spite of the fact that 
the wave functions are different. 

The proton pairing energies of the nuclei ^^^Yb and -^^^Yb show an anomalous 
behavior: For these nuclei, the proton pairing correlations vanish completely 
in the HFBe approach. Taking into account the behavior of the Er isotopes 
in Table 3, it would be expected to obtain for the nuclei ^^^Yb and ^^^Yb 
approximately —4 MeV for the proton pairing energy. In order to understand 
this quenching we have performed a calculation for ^^^Yb similar to the one 
done for the nucleus ^^Zr switching off progressively the Coulomb pairing 
contribution. In Fig. 3a) the proton and neutron pairing energies are displayed 
versus the r] parameter. As it can be seen from this flgure, increasing the rj 
value up to rj — 0.8 produces a linear increase of the proton pairing energy. 
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Fig. 3. a) Variation of the pairing energies versus the rj parameter for the nucleus 
164 )-,^ Proton single particle energies £ (in MeV). The Fermi level is the dashed 
line, c) Occupations of the levels in the canonical basis, versus the rj parameter, for 
the nucleus ^^^Yb. The levels have been labeled by the quantum numbers [nuz Iz]^ 
in an axial basis. 

However, for rj = 0.9 a sharp increase is obtained, going to the unpaired 
solution. By extrapolation of the small rj values to rj = 1.0 the expected 
proton pairing energy for this value would be around —2.5 MeV. However, 
the selfconsistent solution at rj — 1.0 corresponds to a non-paired proton 
system, which seems to indicate that the mean field approximation breaks 
down collapsing to the non-correlated regime. In Fig. 3b) we show the proton 
single particle energies in the canonical basis versus the r] parameter. Here 
we can observe the effect of the sharp phase transition on the single particle 
energies. In Fig. 3c) we can see how the occupations of the levels above (below) 
the Fermi level get empty (full). This sharp phase transition will explain the 
anomalous behavior of these two nuclei. The phase transition itself has more 
to do with the mean field approach than with the exchange terms. This can 
be seen very clearly in Fig.4a) where we show potential energy curves as a 
function of the proton fluctuation for different rj values. That means we have 
solved the HFB equations with an additional constraint on the w.f., namely the 
fluctuation in the proton number. For ^^^Yb and r/ = we find a well defined 
minimum around the superfluid solution, ((AiV)^)^ = 5, i.e., the mean field 
approach makes sense. For rj = 0.9, however, the energy surface is quite flat 
indicating that a theory able to include further correlations would be more 
appropriate. For ^^^Yb, where no pairing collapse takes place, we find well 
defined minima for all rj values, see fig. 4b). 

In order to investigate the high spin properties of these nuclei we have solved 
the cranking HFB equations, eq. (2), in both approximations for the same 
nuclei. In Fig. 5, we display the pairing energies, —E^p, as a function of the 
angular momentum. We find large differences between the approximate and 
the exact results. The absolute values of the HFB^ proton pairing energies 
are much larger than the exact ones, the difference being not only quantita- 
tive but also quahtative for the nuclei ^^^Yb and ^^^Yb. These nuclei remain 
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Fig. 4. Binding energy versus the constrained proton fluctuation for different r] 
values. 



Nucleus (EXP) (HFB^) E^ (HFB,+) E^ (HFBe) A£;^(%) 





91.4 


82.8 


67.7 


68.4 


21.1 


166EJ. 


80.6 


77.6 


63.5 


65.0 


19.4 




79.8 


73.0 


59.4 


60.7 


20.3 


162Yb 


166.9 


115.7 


90.3 


93.9 


23.2 


164Yb 


123.4 


98.7 


85.5 


80.8 


22.2 


166Yb 


102.4 


88.8 


75.9 


75.8 


17.2 



Table 4 

Experimental (1^2 Yb [14], i64Er and ^^^Yb [15], i^^'Er and I'^^Yb [16] and ^^SEr 
[17]) and theoretical transition energies from I = 2h to I = Oh obtained with the 
methods HFB^, the perturbational one HFB^.)., the exact calculation HFBe (in keV) 
and the relative variation between HFB^ and HFBg (%). 

super-fluid for the whole spin range in the approximate solution and normal 
fluid in the exact one. For the neutron system, the behavior is qualitatively 
similar in both approaches, the small differences observed depend on the nu- 
cleus considered. The strong quenching of the pairing energies obtained in 
the exact calculations will affect some observables, for example, the transition 
energies along the yrast band. In Table 4 we show the E^(2+ — > 0+) values 
obtained in the approximations HFB^, HFB^.!. and in the exact calculation. 
As it can be observed, to consider the exchange Coulomb term and all the 
contributions to the pairing field leads to a decrease in the transition energies, 
with a relative change around 20%. Since the perturbational calculation up 
to first order (HFBs+) provides a good approximation to the total binding 
energies we expect that the transition energies calculated in this way will also 
be close to the exact ones. As we can see in the table this is the case. On the 
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other hand, by comparing our results with the experimental ones, which are 
included in the same table, we observe that the HFBe method gives a worse 
agreement with the experiment than the HFB^.!. one. 




12 18 24 
l(R) 



12 18 24 
l(R) 



12 18 24 30 
l(R) 



Fig. 5. El^^[MeV] versus angular momentum calculated with HFB^ (full symbols), 
HFBe (empty symbols). The circles (boxes ) are for protons (neutrons). 

It is interesting to know what happens with the Dl parameterization of the 
Gogny force. We have done the calculation with this parameterization for the 
nucleus ^^''Er, we obtain the value = 95.7 keV with HFB, and = 79.0 
keV with HFBe, with a relative variation between both calculations of 21.1 %, 
the same relative variation as obtained with the DIS parameterization, as we 
can see in the Table 4. The two values of the transition energy E-y obtained with 
the Dl are closer to the experiment, but again the approximate calculation 
HFBg gives the best agreement with the experiment. 



The behavior of the gamma ray energies at high spin is displayed in Fig. 6 for 
the six analyzed nuclei, together with the experimental data. We find that the 
results of the approximation HEB^ are rather different from the exact ones 
as one expects from the disparity of the pairing energies found in the two 
calculations. For the nuclei ^^^Er, ^^^Er and ^*''^Yb we have also included the 
approximation HFBs+. As expected the agreement with the exact calculation 
is much better than the one obtained in the HEB^ calculation. We have also 
checked the Slater approximation for deformed nuclei and at high angular 
momentum. That means, we have evaluated the Eock exchange Coulomb term 
in the Slater approximation and added it to the energy Es of eq. (10) before 
the variation. The results of these calculations, quadrupole moments, pairing 
energies and gamma ray energies along the Yrast band, practically coincide 
with the results of the plain HEB^ calculation. 
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Fig. 6. E^[MeV] versus angular momentum calculated with HFBg (full circles), 
HFBe (empty) and the experimental values (full squares). The empty squares in 
the panels of the nuclei ^^^Er, i^Sgj. ^^'^Yh represent the HFBs+ approach. The 
experimental values are from: ^^^yb [14], ^^-^Er and ^^^Yb [15], ^^^Er and ^^^Yb [16] 
and i^^Er [17]. 

From the results discussed up to now, it seems that it is enough to consider 
the neglected terms in first order perturbation theory. This, however is not 
quite true since some observables can not be expressed as a function of the 
binding energies. For instance, the transition probabilities, magnetic moments, 
nuclear radii, etc. It is interesting to investigate the predictions for some of 
these observables in the exact and in the approximated calculation. In Fig. 7, 
we show the beta deformations as a function of the angular momentum for 
the six nuclei under study. The beta deformation parameter is mainly af- 
fected by the alignment effects and the pairing correlations. In general, larger 
pairing correlations imply smaller deformations and larger alignments smaller 
/3-deformations (Coriolis anti-stretching effect). In the exact calculations, in 
general, we have smaller absolute values of the pairing correlations and larger 
alignment. In principle, one could expect a compensation and only in cases 
where one of the effects is much larger than the other one we should obtain 
a change in deformation corresponding to the largest effect. In the figure we 
observe that in the nuclei ^^^Er, ^^^Er and ^^^Yb such a compensation takes 
place. In ^^^Er the compensation takes place for small spins but not for large 
ones and in ^^"^Yb and ^^^Yb there is no compensation at all. Looking in Fig. 6, 
we find that in ^^^Er, there is clearly a stronger alignment in the exact solu- 
tion than in the approximated one, which causes an smaller deformation in 
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Fig. 7. The /3 deformation parameter versus the angular momentum calculated with 
HFBg (full symbols), HFBg (empty symbols). 

the exact calculation at high spins. In ^^^Yb and ^^^Yb the large quenching 
observed in the pairing energy explains the distinct behavior of the exact so- 
lution from the approximate one. In Fig. 8 we present the reduced transition 



2.4 
2.0 
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12 18 24 
l(R) 

Fig. 8. Reduced transition probabilities versus angular momentum calculated with 
HFBs (full circles), HFBg (empty) and the experimental ones (full squares). The 
experimental values are from ^^^Yb [14], ^^'^Er and ^''^Yb [15], ^^^Er and ^^^Yb [16] 
and i68Er [17]. 

probabilities B(E2) for the nuclei we are analyzing. Though these probabilities 
are related with the /3-deformations we expect larger differences between the 
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two calculations in the B(E2) case because this quantity has to do only with 
protons while to the /^-deformations both the proton and the neutron systems 
contribute. We should remember that neglecting the Coulomb exchange terms 
affects only the proton system. We find this supposition to be right, specially 
at high spins. For the nucleus ^^^Yb there are differences even at zero angular 
momentum due to the change of deformation in the ground state (see Table 3). 
Though we find some differences between both calculation, in the comparison 
with the experimental results it does not matter which one we take. 




Fig. 9. Gyromagnetic factors versus angular momentum calculated with HFEg (cir- 
cles), HFBe (squares). Full (empty) symbols are for protons (neutrons). 

In Fig. 9 we finally present the gyromagnetic factors versus the angular mo- 
mentum. For this observable we expect only small differences between the two 
calculations for the neutron gyromagnetic factor, for the same reasons as 
before, and somewhat larger for the proton gyromagnetic factor gp. We find 
this supposition to be right. In the HFB^ approximation, the general behavior 
of these nuclei is, first, a more or less sharp neutron (ii3/2) alignment take 
place at medium spin values and then a proton (hii/2) alignment take place 
at high spins. In the HFBg calculation the situation is qualitatively similar 
with the exception of the nucleus ^^^Er which shows a very delayed neutron 
alignment, contrary to the experimental situation [16]. 



4-2 Actinide region: The nucleus ^^°Pu. 



In this subsection we shall now discuss an example of a heavy deformed nu- 
cleus, to see if there is a possible mass dependence of the terms neglected 
in the approximate calculations. As a representative of the actinide region 
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we have chosen the nucleus ^ °Pu. In the calculations we have used the basis 
q= 1.3, p= 1.0 and Nq = 14.0 which is large enough to guarantee the conver- 
gence, and the DIS parameterization of the force. Let us first discuss the wave 
function content of the HFB^ approximation and of the exact theory HFBg. 
In Fig. 10, panel (a), we show the pairing energies E^^ in both theories as 
a function of the angular momentum. As before, we find a strong quenching 
for the values of the proton system in the HFBg as compared to the HFB^ 
one. The values for the neutron system, on the contrary, are very similar in 
both approaches. In panel (b) we display the electric quadrupole moment. At 
small angular momentum both predictions are very close but at large angular 
momentum they differ by up to 2% at the largest spin considered. The fact 
that the HFBg prediction is smaller can be understood looking at panel (c) 
where we show the gyromagnetic factors. As expected the gn are very sim- 
ilar in both approaches for all spin values. The gp, in the HFBe prediction 
indicates a larger proton alignment than in the HFE^ one, which will cause 
the observed anti-stretching effect in the HFBg approach. With respect to ob- 




Fig. 10. Properties of the nucleus Pu along the yrast band in the HFB^ and HFBe 
approaches, a) Pairing energies, b) Charge quadrupole moments, c) Gyromagnetic 
factors, d) Gamma ray energies, e) Moments of inertia J'^ and f): Moments of inertia 
J^'^\ The labels for panels a), b) and c) are the ones of panel c). The labels for d), 
e) and f) are listed in d). The experimental data are from [18]. 



servables related to energy differences, we can consider the approaches HFE^, 
HFBs+ and HFBe. In the same figure, we present the transition energies E^, 
panel d), first moment of inertia J'^, panel (e), and second moment of inertia 
JT"^, panel (f), in the three approaches, as well as the experimental data. The 
agreement of the plain HFB^ with the experiment for the three observables is 
spectacular. The exact solution HFBe, according to the less pairing-correlated 
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wave function, see panel (a), provides too large moments of inertia and too 
small gamma-ray energies. For the perturbative approach HFBs+, we obtain 
values similar to the exact calculation HFBe, though for the more sensitive 
quantity J"^ some differences are found at the smallest and the largest angular 
momenta. 

It is very remarkable that the HFBs+ binding- and 7-ray energies are very 
close to the HFBe ones. This fact also happens in the other nuclei calculated 
so far. To gain some insight in the reasons for this agreement we have cal- 
culated separately all contributions to the total energy in both approaches. 
That means, with the wave function determined by minimization of Es-, 
we have evaluated the Fock and pairing term of the Coulomb force, as well 
as the pairing terms stemming from the spin-orbit and the two-body kinetic 
energy terms. We have also separately calculated the mentioned terms with 
the wave function |$)e determined by minimization of the exact energy E^^. In 
the upper part of Table 5 we show the energy Eg of Eq. (10) evaluated with 
the wave function (approach s) and with |$)e (approach e), as well as the 
terms just mentioned (we use the notation of Eq. (9)) for the ground state of 
^^'^Pu. The last column, denoted £'^4, corresponds to the total energy, E^^ in 
the s approach and E^. in the e one. As expected the binding energy is lower 
in the exact than in the perturbative calculation. On the other hand, since the 
wave function |$)s has been determined by minimizing E's, it is also obvious 
that this quantity is lower in the s approximation than in the e one. Concern- 
ing the other terms, each of them separately provides a bigger energy lowering 
in the e than in the s appoach. We also observe that the relative largest dif- 
ference between both approaches corresponds to the pairing Coulomb term. 
In the lower part of the same Table we show, as a function of the angular 
momentum, the energy differences = ^(/) — ^(/ — 2), where ^ stands for 
each energy entry of the upper part of the Table, i.e., represents the con- 
tribution of each term to the gamma ray energy from the state / to 1 — 2. 
These energy differences have been calculated again in the approaches s and 
e. Let us first concentrate on the I = 2h row. The prediction for the gamma 
ray energy taking into account only the terms of the standard approach are 
almost identical independently of the fact that we evaluate them with 
or with |$)e and the same happens with the other terms. Obviously, if one 
adds all these contributions, the prediction for the total 7-ray energy with all 
terms is very similar in both approaches. It is also interesting to note that the 
contributions from the terms not considered in the standard approach arc all 
negative and lead to a more compressed spectrum. We also observe that the 
largest contribution by far is the one from the Coulomb pairing term, though 
the other contributions are not at all negligible. At higher angular momentum 
we observe the same trend as for / = 2h, i.e., the energy differences for a given 
/ are very similar in both approaches. It seems therefore that the reason, for 
the good agreement of the gamma ray energies in both approaches, is that 
the dependence of the exchange terms from the angular momentum is almost 
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Table 5 



Upper part: Ground state contributions to the binding energy of the nucleus ^*°Pu 
in the HFEs^. and in the HFBg approaches. The energies are given in MeV and the 

angular momenta in h. Lower part: Contributions to the gamma ray energies along 
the Yrast band in keV in the HFBs+ and in the HFBg approaches. 

identical in the HFBs+ and in the HFBg approaches. For I = Oh the binding 
energies differ only by about 500 keV. The exchange terms are rather different 
in both approaches, however, when the 7-ray are calculated, the differences 
cancel and one gets similar values in both approaches. 



5 Superdeformed nuclei in the A Ki 190 region 



This region of superdcformation has already been studied using the HFE^ 
method with the Gogny force DIS [6,7]. In our calculations we use the basis 
q= 1.5, p= 1.0 and Nq = 12.5 and the DIS interaction. Taking into account 
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the results obtained in the other regions, we do not expect a big influence of 
the Coulomb exchange terms because for the superdeformed states the proton 
pairing energies either vanish or are very small even with the HFB^ method. 
A large change of the pairing energies due to the Coulomb exchange and the 
other terms cannot exist, therefore, in this states. As an example of superde- 
formation we have chosen the nucleus ^'^'^Hg. From this example we will learn 
how the other pairing contributions (spin-orbit and two-body correction) af- 
fect the neutron pairing energies and the other nuclear properties through the 
selfconsistency. As before, we shall first look for the wave function content of 
the HFBs and the HFBe approaches. In Fig. 11, panel (a), we can see the 
behavior of the pairing energies in both methods. As mentioned, the proton 
pairing energies along the band vanish in both calculations. The neutron pair- 
ing energies are close to each other in both calculations, though somewhat 
smaller values are obtained for the HFBg at low angular momentum. The 
charge quadrupole moment is depicted in panel (b). Only small differences 
between both approaches are found at large angular momentum. The gyro- 
magnetic factors, panel (c), are again very close to each other. We shall now 
turn to a comparison between the HFB^, HFBs_|_ and HFBe predictions for 
the 7-ray energies, panel (d), first moments of inertia, panel (e) and second 
moments of inertia, panel (f). For these observables we do not find remarkable 
differences between the three approaches. From this section we may conclude 




that the Coulomb exchange and pairing contribution barely modify the prop- 
erties of the nucleus when the approximate calculation HFE^ has vanishing 
proton pairing energies. Then, the approach RFB^ is a good approximation to 
the completely self-consistent calculation, HFBe, even for the wave function. 
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6 Energy Surfaces 



In this section we shall finally investigate the effect of the exchange terms 
as a function of a collective variable at zero angular momentum. We have 
used an axially symmetric HFB code that allows reflexion asymmetric shapes. 
We have performed two different calculation. In the flrst one we have ne- 
glected all exchange terms, i.e., the HFBs approach, but we have included, 
the Fock Coulomb term in the Slater approximation, we shall call this ap- 
proach HFBs,s/a- The second approach includes all exchange terms with the 
exception of the contribution of the spin-orbit term to the pairing fleld,|3] we 
shall call this approach HFBe^nso- Notice that since we do not add at the end 
of the calculation HFBg^s/a the pairing terms of the Coulomb and the two- 
body kinetic energy (both contributions are repulsive) it is possible to find 
HFBs solutions with an energy deeper than the HFBg approach. That 
means we are more interested in the general behavior of the energy surface 
than in the absolute values. The most common case of energy dependence as 
a function of a collective variable appears in the calculation of fission barriers. 
In this calculations we use the basis q= 1.5, p= 1.0 and Nq = 15.1 and the 
DIS parameterization. In fig. 12, on panel a), we display the binding energy of 
^^^No as a function of the constrained quadrupole moment in both approaches. 
We find a similar behavior in both calculations, though small differences are 
noticeable. For example, the height of the first barrier is slightly larger in the 
IIFBe,nso approach than in the HFBs^s/a. The second one, however, has more 
or less the same height in both approaches. In panel b) we display the pair- 
ing energies along the fission path. The neutron pairing energy is practically 
the same in both calculations. The proton pairing energies in the HFBg nso 
calculations show the Coulomb antipairing effect already discussed in other 
nuclei. They are in absolute value 5 to 10 MeV smaller than in the calculation 
without Coulomb pairing term. 

We shall now investigate the effect of the exchange terms in cases where the 
reflexion symmetry is broken i.e. the wave functions do not have a good parity. 
In these calculations we use the same theoretical approaches as in the ^^^No 
ones, we use the basis q= 1.3, p= 1.0 and Nq = 13.1 and the DIS interaction. 
In panel c) of the same flgure, we present the binding energy of the nucleus 
^^^Th versus the constrained octupole moment. As with the quadrupole oper- 
ator we do not flnd any remarkable difference between both calculations. The 
depth and shape of the octupole deformed minimum are almost identical in 
both approaches. In panel d), flnally, we present the dipole moment versus 
the constrained octupole momentum. Again, small differences between both 



We do not include this term in the calculations just because it is not relevant 
for the problem we are considering, besides the fact that it is very cumbersome to 
include in the axial HFB code. 
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Fig. 12. a) Results in the HFB^ approach, crosses, and HFBe,nso) filled circles, 
a) Binding energy of the nucleus ^■'^'^No versus the quadrupole moment, b) Pairing 
energies along the fission barrier, protons (neutrons), dashed (solid) line, c) Binding 
energy of the nucleus 226Th versus the octupole moment, d) Dipole moment versus 
the octupole moment. 

calculations are found but without further relevance. 



7 Discussion 

The main outcome of our investigation is that the exchange terms usually 
neglected affect mainly the nuclear pairing properties. More specifically, they 
affect mostly the proton nuclear pairing properties because of all missing terms 
the largest contribution to the pairing field arises from the Coulomb potential. 
Observables little or not dependent on the pairing correlations, are, therefore, 
not affected by these terms. On the contrary, observables like the moments 
of inertia or two-quasiparticle states energies, strongly dependent on pairing 
correlations, are very sensitive to these terms. 



We have found that the results of calculations performed with the plain HFB^ 
approximation differ significantly, for some observables, from those of the exact 
solution HFBg, in nuclei where the proton pairing energies are relatively large. 
On the contrary, the HFB5+ results agree, in general, very well with those of 
the exact solution. In cases of very small proton correlation energies all three 
approaches, HFBg, HFB^.!. and HFBg practically coincide. 

For those situations where proton correlations are relevant, the predictions of 
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the HFBs approximation agree better with the experimental data than the 
exact ones or than the ones from the HFBs+ approximation. Paradoxically, 
the Gogny force was fitted [2] using the exact HFB approach to calculate 
properties of spherical nuclei. The pairing properties of the Gogny force were 
adjusted by fitting the odd-even mass difference of the Sn isotones. The Sn 
isotones, however, have a major shell closure in protons, Z— 50, and in this 
case, practically, it docs not matter whether one performs approximate or 
exact calculations, i.e., effectively the Coulomb exchange terms were not con- 
sidered in this fit. Furthermore, the fact that including the exchange terms 
in the calculations worsens the good agreement of the theoretical results with 
the experimental data, reinforces the conjecture that the actual fits of the 
Gogny force should be used, in calculations with open proton shells, without 
the mentioned terms. 

Of course the Gogny force is isospin invariant but the Coulomb force is not. 
The Hartree-Fock potential that one obtains after solving the HFB equations 
neither is isospin invariant. It is not clear, therefore, that a different parame- 
terization (concerning the pairing part of the interaction) of the force would 
not be obtained if nuclei with open proton shells had been considered in the 
fit. The whole issue of adjusting the pairing properties of any force is itself 
rather fuzzy. Most of the fits have been done to reproduce the experimentally 
observed odd-even mass difference. This energy difference can, however, not 
be attributed only to pairing correlations [20]. From the theoretical point of 
view the calculation of this energy difference presents also problems because 
binding energies of odd nuclei are not easy to calculate. Blocking effects and 
angular momentum projection should be performed in order to evaluate prop- 
erly this energy, besides additional couplings. However, in order to save CPU 
time during the fit, approximations are used, in particular in order not to 
break the time reversal symmetry the equal filling approximation is used. 

Taking into account the complication of the neglected exchange terms, see 
Appendix B, the fact that the HFBg approximation, commonly used in most 
HFB calculations, provides such good agreement with the experimental data 
is very fortunate. 



8 Conclusions 

For the first time we have performed an exact self-consistent HFB calculation 
in a triaxial basis with the Gogny force. The exact Coulomb exchange term and 
the contributions of the spin-orbit, two-body correction of the kinetic energy 
and Coulomb terms to the pairing field have been included in the calculations. 
An exhaustive analysis of each term has been performed for spherical nuclei 
in different mass regions. We find the Coulomb contribution to the pairing 
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field to be the most relevant one. The other terms, though not negligible, give 
smaller contributions. 

For deformed nuclei we have studied the commonly used HFB^ approximation, 
the perturbative HFBs+, and the exact sohition. We have been concerned with 
different mass regions at zero and high spins. As in the case of spherical nuclei, 
we observe a reduction of the proton pairing energies in the exact calculations 
as compared to the HFBs ones. Ground state energies are very similar in 
the perturbative and in the exact calculations. Energies of excited states are 
somewhat different in the HFB^ approximation than in the exact one. We have 
also investigated superdeformed states. Since the proton pairing correlations 
of these nuclei arc small we do not find any significant difference between the 
approximate and the exact solution. Lastly, we have analyzed the effect of 
the most relevant exchange terms in the calculation of energy surfaces. The 
energy as a function of the quadrupole (octupole) mass operator is practically 
identical in the approximate and in the exact calculations. 

We also have found, that the Slater approximation for the Fock Coulomb term 
is a good one for all HFB calculations, in agreement with ref.[3] for the HF 
case. 

For all nuclei and situations analyzed the following general conclusions apply: 

(1) In cases where the proton pairing correlations do not play an important 
role, anyone of the approximations considered practically coincide with 
the exact calculation. 

(2) When proton pairing correlations are relevant, then: 

• For energy-related observables the perturbative approach HFB^.!. is a 
good approximation to the exact one. 

• For those observables whose evaluation requires the wave functions, one 
should perform the exact calculation. 

• If the HFB wave functions are thought as a basis for more elaborated 
theories, like RPA or GCM, then the exact calculation is required. 

(3) The way in which the fits of the Gogny force were performed favors 
neglecting the mentioned exchange terms in HFB calculations using these 
fitted parameter sets. 

It remains to investigate the real importance of this terms by performing new 
fits of the Gogny force taking into account explicitly all terms discussed here. 

We thank J.F. Berger, M. Girod and S. Peru for discussions. This work was 
supported in part by DGICyT, Spain under Project PB97-0023. One of us 
(M.A) acknowledges a grant from the Spanish Ministery of Education under 
Project PB94-0164. 
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Appendix A: The interaction 



As an effective interaction we use tlie Gogny force [2] 



2 



vi2 = E e-(^"*^-^"*^)'/'^'(W, + BiP, - HiP, - MiP^P,) + 




(-1) 



and tlie Coulomb force 




(l + 2ri,)(l + 2r2,) — 



n - r2 




(-2) 



Additional contributions taken into account in the calculations arise from the 
one-body and two-body center of mass corrections 



In the calculations we use the parametrizations DIS [8] and Dl [2,9]. 

Appendix B: Calculation of the Coulomb Hartree-Fock field and pairing tensor. 

To compute the Hartree-Fock field and pairing tensor for the Coulomb inter- 
action we have followed the idea used in [4]. It amounts to use the identity 



to compute the Coulomb matrix elements in terms of those of Gaussians with 
range /i. The main advantage of this approach is that the matrix elements of 
the Gaussians are factorizable in terms of quantities defined for each Cartesian 
direction x, y and z. Unfortunately, the price to pay is the // integration. 
Introducing the HF field ^kik2il^) pairing tensor A'^^^^^i), see eq.(6), 

of a Gaussian of range fj, we can express the corresponding quantities for the 
Coulomb potential as 




(.3) 
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and 

To transform the integration interval to a finite one, one makes the change of 
variables = (1 + a//^) ^, with a a parameter to be determined below. We 
finally obtain 

pC _ 0,2 / du a ( ^^{l-u^) \ 

^ '^'^ ~ y TT I (1 - U^f [ U ) 

and 

The integrations arc carried out numerically with the Gauss- Legendre method. 
The optimal choice of a, turns out to be equal to the maximum of (a^, Oiy, a^) 
where 0^1 = ^ and 6j are the oscillator lengths (see ref. [4]). 

i 

In the numerical computation of the Coulomb exchange and pairing terms we 
have used 20 points for the Gauss-Legendre integration. This number of points 
has proven to be enough for a precise determination of the exchange field and 
pairing tensor. The number of integration points used makes the evaluation of 
the Coulomb terms rather costly in terms of CPU time: it usually takes a factor 
5.5 longer to compute those factors than to calculate the Fock and pairing parts 
of the Brink-Boecker interaction (in this case, we have two gaussians and we 
have to compute them for both protons and neutrons). The calculation of the 
exchange and pairing parts of the BE interaction is already the most costly 
part of each iteration as it usually takes between a factor 10 and 15 longer 
than the evaluation of the direct terms. In a standard workstation and for 
a configuration space with parameters p = 1, q = 1.3 and Nq = 11.1 the 
evaluation of the direct field takes 1.9 seconds, the combined Fock and pairing 
terms of the BB part takes 19 seconds and the evaluation of the Fock and 
pairing parts of the Coulomb interaction takes 103 seconds. One has to keep 
in mind that the previous numbers are for each iteration of the minimization 
process. If we increase the size of the basis by using g = 1.6 and A^o — 14.0 the 
previous numbers climb up to 7, 92 and 540 seconds respectively. We clearly 
see the tremendous slow down introduced in the calculations stemming from 
the consideration of the Fock and pairing terms of the Coulomb interaction. 
This tremendous slow down can be even worse for the practitioners of the 
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Skyrme interaction as their CPU time consumption per iteration has to be of 
the order of the one needed to calculate the direct field of the Gogny force. 

Concerning the computation of the pairing field associated to the spin-orbit 
interaction we have followed the ideas aheady presented in [4]. An explicit 
expression for all these terms can be found in [21]. The numerical computation 
of this field takes almost the same time as the computation of the HF field 
and therefore, its inclusion is not significant in terms of CPU time. 
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